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Abstract 

The observed branching ratios for B Kr]' decays are much larger than 
factorization predictions in the Standard Model (SM). Many proposals have 
been made to reconcile the data and theoretical predictions. In this paper 
we study these decays within the SM using flavor U(3) symmetry. If small 
annihilation amplitudes are neglected, one needs 11 hadronic parameters to 
describe B PP decays where P can be one of the n, K, rj and ry' nonet 
mesons. We find that existing data are consistent with SM with flavor U(3) 
symmetry. We also predict several measurable branching ratios and CP asym- 
metries for B K{iT)ri(rj'), 'r]{i]')r]{ri') decays. Near future experiments can 
provide important tests for the Standard Model with flavor U(3) symmetry. 
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Experimental data from CLEO, BaBar and Belle [1-4] have measured branching ratios 
of S — > Kri' around 6 x 10~^ which are substantially larger than theoretical calculations 
based on naive factorization approximation in the Standard Model (SM) [5]. Although there 
are some improvements in calculating the branching ratios in the last few years by using 
QCD improved factorization method [6] , there are still large uncertainties in calculating the 
branching ratios for B — > Kr]' because of issues related to 771 — r]^ mixing and QCD anomaly 
associated with 771. There are also many speculations about possible new physics beyond the 
SM in these decays [7]. Before any claim can be made about new physics, one must study 
the SM contributions in all possible ways to see if it is really inconsistent with experimental 
data. 

In this paper we carry out a systematic study of S — > Kr]' , and more generally processes 
of S — > PP decays with at least one of the P to be 77(77') in the final states by using 
flavor symmetry in the SM. This way one can relate different decays to predict unmeasured 
branching ratios and CP asymmetries. Drastic deviations between the predicted relations 
and experimental data can provide information about the SM and models beyond. Similar 
considerations based on SU{?>) have been applied to S — > PP decays [8], and shown to be 
consistent with data [9]. If one sticks to flavor SU(3) symmetry, one needs to introduce the 
singlet 771 in the theory and to add additional amplitudes to describe new decay modes [10]. 
One may also consider to promote the flavor SU(3) symmetry to flavor U(3) symmetry such 
that 771 is automatically included in the theory. 

Flavor U (3) symmetry has been studied in Kaon decays. There there are non-negligible 
symmetry breaking effects. For B decays one may also expect symmetry breaking effects 
to exist. There are also some studies of U(3) symmetry for B decays [11]. Present data, 
however, are not able to make clear statement about whether this symmetry is badly broken. 
In this paper we will take the flavor U(3) symmetry as working hypothesis and to study 
whether experimental data can be explained by carrying out a systematic analysis. We flnd 
that the SM with flavor U(3) symmetry can explain all existing data, in particular can obtain 
large branching ratios for B Kr)' decays. We also predict some unmeasured branching 
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ratios and CP asymmetries which can be used to further test the theory. 
The quark level effective Hamiltonian can be written as [12] 
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Hlff = ^[K.v;;(ciOi + C2O2) - E(K6KVr + Vt,v,\cf)oi (i) 

Here Vij are KM matrix elements. The coefficients Ci_2 and c?^ are the Wilson Coefficients 
which have been evaluated by several groups [12] with |ci,2| >> |cj''|. Oi are operators 
consist of quarks and gluons. 

The B — > PP decay amplitudes can be parameterized as 



(2) 



A{B ^ PP) =< PP\HyB >= -^[VubV:^T + Vtt,V*P], 

where B — {Bu,Bd,Bs) — {B^,B°,Bg) which form a fundamental representation of 
SU{3) (or {U{3)). The amplitudes T and P are related to the hadronic matrix elements 
< PP\Oi\B > which are very difficult to calculate. For our purpose, however, we only 
need to know the fact that under SU{3) (or C/(3)) Oi,2, 03_6,ii, and O7-10 transform as 
3 + 3' + 6 + 15, 3, and 3 + 3' + 6 + 15, respectively [8], and to parameterize the amplitudes 
according to SU(3) (or U(3)) invariant amplitudes to be discussed below. 

As mentioned earlier that there are two approaches to the problem from flavor symmetry 
point of view. We first work with the U(3) symmetry approach. In this case, the tt, K, rjs 
and rji form a nonet representation of U{3), Mj, with 
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One can write the T amplitude for B 
amplitudes as 



TS^i 
PP in terms of the flavor U(3) invariant 



V 



T = AlBiH{3)\M^M[) + CiBiMlM}H{3y 
+ AlB,H{QriM]M^ + ClB,M]H{QtMi 
+ A^B.HiW^MW + C^B,Mm{T5)j'Mi 



+ B^BiH{?,YMjMl + BlBiH{<ofiM^Ml 

+ B^BiH{^yiM^M\ + DlBiM]H{?>yMl . (3) 

In Table I we list all decay amplitudes involving 771 g. The amplitudes containing only K 
and TT in the final states can be found in Ref. [9]. There are a few new features for the U(3) 
amphtudes in eq. (3) compared with the SU(3) amplitudes for B — > PP. The last four terms 
are new. In SU(3) case because the traceless condition of M/, these terms are automatically 
zero. With SU(3) symmetry, the amplitudes Aq and Cq always appear in the combination 
of — Aq [8] . This degeneracy is, naively hfted in processes with 771 in the final states. It 
seems that there is the need of having both Cq and Aq to describe the decays increasing the 
total number of hadronic parameters by one. However, this is not true since that the Aq 
and Cj terms in decay modes with rji in the final state can be written as = Cj — Aq, 
and the additional Aj be absorbed by redefining the amplitude Bj — Bj + Aj. In Table 
I we therefore have listed the decay amplitudes in terms of the independent U(3) invariant 
amphtudes, Cj^^^, AT-, BT^-^ and Dj. 

We now describe the other approach to include r^i in S — > PP decays. Here one treats rji 
as a singlet of SU (3) and parameterizes the decay amplitudes according to SU(3) symmetry. 
In this case there are also additional four new terms, 

Tnew = a^BiH{3y7]i7i, + b^B,M;{8)H{3y7i, 

+ c^BiH{<o)fMi{%)r], + (fBiH{TE)fMl{8)ri,. (4) 

In the U(3) symmetry limit, we have 

= Ai + 353^ + ^CT + Dl IF = ^CJ + VSDl 

We also see from the above that one can use Cj and Bj to absorb Aq by writing, = 
V^Bq + Cj. It is interesting to note that both approaches discussed above introduce the 
same number of new parameters, four of them, into the theory. In our analysis we will work 
with the fiavor U(3) symmetry. 



To obtain the amphtudes for B decays with at least one 77(77') in the final states, one also 
needs to consider r] — rj' mixing, 
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The average value of the mixing angle 9 is —15.5° ± 1.3° [13]. We will use 9 — —15.5° in 
our fit. 

There are similar U{3) invariant amplitudes for the penguin contributions. We indicate 
them as te-, ^Cte-, and D^^ The amplitudes Ai and Bi are referred as annihila- 

3,6,15' 3,15' 3,6,15 3,6,15 ^ ' ' 

tion amplitudes because the B mesons are first annihilated by the interaction Hamiltonian 
and two light mesons are then created. In total there are 22 complex hadronic parame- 
ters (44 real parameters with one of them to be an overall unphysical phase). However 
simplification can be made because the following relations in the SM, 



^6^(^6^) - ~ 
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cj(sj) ^ -omsc^iB^) , 

-CU^h^ B^) ^ +0.015C§(^^, Ei^). (7) 



2 ci + C2 - 3(cr + cYS)/2 
Here we have used the Wilson Coefficients obtained in Ref. [12]. 

We comment that in finite order perturbative calculations the above relations are renor- 
malization scheme and scale dependent. One should use a renormalization scheme consis- 
tently. We have checked with different renormalization schemes and find that numerically 
the changes are less than 15% for different schemes. In obtaining the above relations, we 
have also neglected small contributions from which cause less than 1% deviations. 

Using the above relations the number of independent hadronic parameters are reduced 
which we choose to be, C^'^ {A^'^)^ Cj, Cjg{Aj^), B^'^, B'^, Bj^, D^'^ . An overall phase 
can be removed without loss of generality, we will set to be real. There are in fact only 
25 real independent parameters for B PP in the SM with ffavor U(3) symmetry, 

C^, CTe^'^ Cje^'^, C^^'T,, ATe'^^ , A^e'^i , A^^^t,, 



15 
i<5. 



iSr- 



B-U 3 , B^e s , BU 6 , Bj^ DU 3, D^e 3 



TABLE I. U(3) decay amplitudes for B PP with at least one of the P being a r/s or rji. 

AS = AS = -1 

T^M) = 75 + + + T-Us) = -^,{C- + CJ - AI^-5C^), 
TnX id) = ^(2Af + iCj - C- -AI^+ C^, T^, (s) = ^(C^ + 2AZ_ - 2C§), 

^Kijd) = -^(^r + CI -A^- 5Cg, T^is) = V2{AT + |Cj -A^- 2C^) , 

(d) = ^si^CT + CI + 6^^ + Tl-_^^ {s) = ^^{2CT + Cj + 6A^ + 

+3BJ + 9Bl^ + +3BJ + bI^ + 3DJ), 

TXM) = ^eC'CT + d - 10AI^-5C^ m^S^) = ^(2C7J - - 2A2_- C7§ 

+3^6^ - 15B§ + 3DT), -3BJ - 3B^ + 3DT), 

id) = ^(2CJ - 3Ci + 6AI^ + 3C§ T^, (s) = ^(C^ - 4Al^ - 2C^ 

-9BI + 9B^ + 3DT), -3Bj + QB^), 

T,%id) = f (3^J + + 9BI + 3D^), T^is) = -#(2Cj - 6A^- 3C§ 
^xo,. id) = ^(2CJ - Ci ~ 2AI^ - -9BI^ + 3Dl), 

-3Bi-3Bl^ + 3DT), T,^^^(a) = f (3AT + + 9Bj + 3Z?J), 



Further the amplitudes and Bi correspond to annihilation contributions and are expected 
to be small which is also supported by data [9]. If the annihilation amplitudes are neglected, 
there are only 11 independent hadronic parameters 

Ci, Cie^'-^ Ce^e*, ^^6^%, D^e^'^r, i^feS^ (8) 

TP T P 

The phases in the above can be defined in such a way that all C/'" and £>■ are real positive 
numbers. 

At present many B — > PP decay modes have been measured at B-factories [2-4]. It 
is tempting to use experimental data to fix all the hadronic parameters described earlier. 
It has been shown that if processes involving r]{r]') are not included, it is indeed possible 
to determine all the SU(3) invariant amplitudes, Ai and Ci [8]. When processes involving 
77(77') are also included, a meaningful determination of all hadronic parameters (25 of them) 
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is, however, not possible at present because of too many parameters. We therefore in the 
following neglect the annihilation amplitudes, which are anticipated to be small, to see if all 
data can be reasonably explained, in particular to see if large B — > Kr]' branching ratios can 
be obtained, with only 11 parameters given in eq. (8). This is a nontrivial task. Remarkably 
we find that all data can, indeed, be well explained. 

We use the averaged CLEO, BaBar and Belle data [2-4] shown in Table III and IV to 
fix the unknown 11 hadronic parameters by carrying out a global analysis. The results 
are shown in Table II. In our analysis due to the lack of knowledge of the error correlations 
from experiments, in obtaining the averaged error bars, we have, for simplicity, taken them 
to be uncorrelated and assumed to obey Gaussian distribution taking the larger one between 
(T+ and (T_ to be on the conservative side. Experimental data on e;^, B — B mixing, |V^b|, 
iKtft/Kbl) and sin 2/5 provide very stringent constraints on the KM matrix elements involved 
in our analysis [9,14,15]. We have treated them in our analysis as known parameters with 
the values A = 0.2196, A = 0.854, p = 0.25 and 77 = 0.34 (7 = 53.4°) determined from the 
most recent data. 

Using the above determined hadronic parameters, we study several other unmeasured 
branching ratios and CP violating rate asymmetries Aqp for B — > PP defined by, 

T{B, ^ PP) - V{B, ^ PP) 



i-CP 



(9) 



r{Bi ^ PP) + r{Bi ^ PP) ■ 

The results are shown in Tables III and IV. 

We now discuss some implications of the results obtained and draw conclusions. We see 
from Tables III and IV that the best fit values for the known branching ratios are in good 
agreements with data, in particular large B — > Kr]' can be obtained. The minimal in our 
fit is 16 with total of 23 data points from B PP decays and the 11 hadronic parameters 
in eq.(8) as fitting parameters. The value of 1.33 for the per degree of freedom represents 
a reasonable fit. In our fit the rj — rj' mixing parameter 9 is fixed at the average value 
determined from other data [13]. We checked the sensitivity of the final results on 6 within 
the allowed region and find the changes are small. We also find that if one reduces the U(3) 



TABLE II. The best fit values and their 68% C.L. ranges for the hadronic parameters. 





central value 


error range 
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"3 


85.6° 


29.8° 




79.0° 


17.4° 




8 9° 


15 3° 




0.122 


0.011 




0.940 


0.340 




-85.0° 


6.0° 




-83.7° 


16.5° 



symmetry to the SU(3) case, just fitting data on S — > tttt, ttK, KK, the values obtained for 
Cj are not very much different than what obtained here. This indicates that the parameters 
Ci are stable when promoting SU{3) to U{3). The large branching ratios for B Kt]' are 
due to the new parameters D^. The U (3) assumption can provide a good approximation for 
B — > PP decays. The Standard Model with flavor U(3) symmetry is not in conflict with 
existing data. 

In our fit, we did not include the branching ratios which only have information on their 
upper bounds, such as Br{Bd K-K+), Br{Bd K^K^) and Sr(S„ n'r]'). Since 
we neglected annihilation contributions, the mode Bd — > K~ has vanishing branching 
ratio which is consistent with data. The stringent upper limit [2-4] of 10~^ on B^ K~K^ 
branching ratio supports the expectation that annihilation contributions are small. The 
predicted branching ratio of 0.7 x 10~^ for B^, K^K^ is safely below the experimental 
bound. The predicted branching ratio for Br{Bu '^~v') is 16.81^1^ X 10"^ The central 
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TABLE III. The central values and 68% C.L. allowed ranges for branching ratios (in units of 
10~^) and CP asymmetries for processes with no rj or rj' in the final states. 





Branching Ratios 


CP Asymmetries 




Exi")(M'i)lK')li 


Fit 


Ex|")('rim(nil 


Fit 


iJu — *■ TT TT 


5.6 ±0.9 


c c+0.9 
'^■^-0.9 


0.06 ±0.16 


0.00 


±Ju — s- A A 


-0.6 ±0.8 






-0.49t0;l^ 


D . _-+_ — 

13(1 — ^ TT TT 


4.8 ±0.5 




0.51 ±0.19 


0.45^0;}^ 


R ^ ^0^0 
-Drf ^ TT TT 


2.0 ±0.8 


1 q+0.8 






±Jd — > A A 




'-'■'-0.2 




-0.49i°;l^ 


— TT A 


18.2 ± 1.7 


2o.illi 


0.04 ± 0.08 


r, r)9+0.03 


-Du ^ TT A 


12.9 ± 1.2 


10.81?:^ 


-0.10 ±0.08 






lo.O □= l.U 


18.9_o4 


n no -I- n n/i 
— U.uy ± U.U4 


n 1 1 +003 

-0.11_o.o3 


Ba ^ ttO^o 


10.3 ± 1.5 




0.03 ± 0.37 


-0 06+° °^ 




















-0.33 


Bs K+K- 




17.8lgi 




_n J 1+0.03 
'-'•^^-0.03 


Bs ^ K^K^ 








^•U^-0.04 



value is slightly larger than the 90% C.L. allowed upper bound 12 x 10"^. But the 68% 
C.L. range is consistent with data. At present it is too early to claim conflict of theory with 
data. But this mode can be used to test the theory. Should a branching ratio much smaller 
than the central value predicted here be measured in the future, it is an indication that the 
assumptions made need to be modified. 

Using existing experimental data, we have determined 11 hadronic parameters needed 
to describe B — > PP decays with flavor U(3) symmetry. We have compared with QCD 
improved factorization calculations and found that the magnitudes of the parameters are 
of the same order of magnitude. In our fit, we determined two U(3) invariant amplitudes. 
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TABLE IV. The central values and their 68% C.L. allowed ranges for branching ratios (in units 
of 10"®) and CP asymmetries with at least one of the final mesons to be a 77 or r( . 





Branching Ratios 


CP Asymmetries 




Experiment 


Fit 


Experiment 


Fit 


Bu^n r] 


4.1 ±0.9 


4 1+0.9 


-0.51 ± 0.20 


-0.23;°;}^ 


Bu ^ K-ri 


3.1 ±0.7 


q q+0.6 

o.o_0.4 


-0.32 ±0.22 


n 07+0.08 


Bu ^ K-7]' 


77.6 ±4.8 


72.8l|9 


0.04 ± 0.04 


07+0-04 

^•^'-0.04 


Bd ^ K% 


2.6 ±0.9 






-0 21+0-07 


Bd ^ 


58.3 ±6.0 


66.511:^ 


0.08 ±0.16 


12+004 


Bu 71"- r/' 
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-U.i»_o.o9 


Bd ^ 7rO?7 




1 9+0.6 




-0.94l°:Ji 


Bd ^ TT^' 








-0.38l°:i 


Bd r)T] 
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"^•-•■-l.l 




_n 00+0.10 
"■''^-O.IS 


Bd r]r]' 




'■0-3.4 




-0 20+0-12 


Bd ^ f?'r/' 








•J- -'--0.28 


Bs^Kt] 




2.8ii:^ 




17+0-07 
^•-'-'-0.07 


Bs ^ Kri' 




19 1+^° 
^^■^-9.0 




-0 S9+0-17 


Bs ^ Tr% 




nf^+o io 

^■^'^-0.10 




0.981H1 


Bs Tr%' 




o.o7l°:i° 




n 01 +0.09 

^•^-'--0.10 


Bs -nr] 




7.ol}:i 




n 90+0.09 
"•^^-0.09 


Bs ^ m' 




24.111J 




0.081°:°! 


Bs ^ 




68.3lt^ 




o.o9l°:°i 
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TP c 

D^' . These are particularly difficult to estimate from theoretical calculations because 
these amphtudes may be related to QCD anomalies. Also in factorization approach, it is 
not possible to reliably calculate the phases in the hadronic parameters. In our fit, we find 
that these phases can be sizeable. Further improved theoretical method is needed to have a 
better understanding oi B ^ PP decays. 

Using the hadronic parameters determined from existing data, we have predicted several 
unmeasured branching ratios. These predictions can be used to test the theory. There are 
six modes involving at least one rj (or a rj') in the final states for Bd decays. Among them 
Bd —>■ 7r~rj' has the largest branching ratio, it is a clear test for the theory. The other 
modes are also in the reach of near future data from B-factories. There are also seven Bg 
decay modes with at least one t] or rj' in the final states. Several of the branching ratios 
are predicted to be large, in particular the predicted branching ratio for Bs —>■ rj'r)' is about 
7 X 10~^ which can be measured at future hadron colliders and can provide another crucial 
test for theory. 

We have also obtained interesting predictions for CP asymmetries in B ^ PP decay 
modes. Many of the predicted central values for the CP asymmetries are larger than 10% 
which can be measured in the near future. These modes can provide important information 
about CP violation in the Standard Model. 

In conclusion, we have carried out an systematic analysis for B — > PP decays in the SM 
with fiavor U(3) symmetry. This approach allows one to study B decays involving at least 
one T) or rj' in the final states. We find that all existing data can be explained, in particular 
large branching ratios for B — > Krj' are possible. There is no confiict between the Standard 
Model and present experimental data. We have also predicted several unmeasured branching 
ratios and CP asymmetries within the reach of near future B-factories. Future experimental 
data will provide crucial information on fiavor symmetries and also the Standard Model. 
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